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ABSTRACT 


The  design,  construction  and  implementation  of  a 
computer  controlled  radio  frequency  (RF)  noise 
measurement  and  recording  system  in  the  30.45  MHz  and 
150.75  MHz  range  and  the  interfacing  of  this  system  to  a 
digital  satellite  link  is  presented.  Earthquake 
prediction  by  the  use  of  various  precursors  and  the 
specific  use  of  electromagnetic  emissions  in  the  RF  range 
as  a  precursor  to  future  earthquake  activity  is 
described. 
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I.  INTRODUCTION 


A.  THE  NEED  FOR  EARTHQUAKE  PREDICTION  RESEARCH 

Intensive  academic  and  social  attention  ha s  been 
drawn  to  earthquake  prediction,  especially  in  countries 
like  Japan,  the  United  States,  China  and  the  Soviet  Union 
which  have  all  suffered  from  very  destructive  earthquakes 
on  numerous  occasions.  Ihvestigation  into  the  realm  of 
earthquake  prediction  in  these  countries  has  brought  out 
clues  pointing  to  the  fact  that  the  prediction  of  some 
earthquakes,  if  not  all  of  them,  will  be  attainable  in 
the  foreseeable  future.  Of  course,  all  of  the  permanent 
measures  that  are  used  to  reduce  earthquake  hazards,  such 
as  better  engineering  design,  improved  land  use  planning, 
and  upgraded  building  codes  and  construction  practices 
will  continue  to  be  needed  even  though  earthquake 
predictions  may  eventually  become  completely  reliable. 
The  benefits  to  society  from  accurate  earthquake 
prediction  will  be  substantial  in  the  savings  of  lives 
and  property. 

It  has  been  estimated  that  the  knowledge  of  an 
impending  great  earthquake  a  year  or  more  in  advance 
would  result  in  a  large  reduction  in  losses.  The  savings 
would  result  from  measures  taken  to  strengthen  buildings 
and  their  contents,  reduce  fire  hazards,  increase  dan; 
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safety,  enhance  nuclear  power  plant  safety  and  so  on. 

For  shorter  term  predictions  such  as  a  week  in  advance, 
substantial  savings  of  lives  would  result  from  temporary 
measures  such  as  the  evacuation  of  dangerous  buildings 
and  the  mobilization  of  emergency  forces.  [Ref.  l:pp. 
35-*i53 

It  is  known  in  many  instances  that  the  earthquake 
itself  will  cause  less  damage  and  casualty  than  the 
events  which  occur  following  the  actual  geophysical 
phenomena.  Fires  and  explosions  resulting  from  ruptured 
gas  lines,  gasoline  storage  tanks,  power  lines  and 
transformers  will  contribute  to  the  severe  loss  of  life 
and  property.  If  the  utilities  such  as  gas  and 
electricity  were  to  be  equipped  with  sensing  devices 
which  would  automatically  interrupt  services  when 
precursory  information  was  acquired  concerning  a 
potential  earthquake,  great  devastation  could  be  spared. 
If  the  ability  to  predict  an  earthquake  could  be  soundly 
established,  a  more  calculated  reaction  to  the  event 
could  be  made.  Evacuation  of  endangered  areas,  or  even 
the  planned  demolition  in  advance  of  structures  deemed 
unsafe  in  the  event  of  a  major  temblor  could  be 
accomplished.  Unfortunately,  the  ability  to  accurately 
and  reliably  predict  the  location  and  intensity  of  a 
forthcoming  earthquake  is  unavailable  at  the  presei.t 
time. 
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Great  excitement  currently  exists  among  seismologists 
over  major  achievements  in  the  effort  to  predict 
earthquakes.  Various  physical  phenomena  precursory  to 
earthquakes  have  been  reported.  Measurable  physical 
precursors  are  now  coming  under  intense  scientific 
scrutiny  and  many  appear  to  be  valid.  Several  scientific 
predictions  have  already  been  successful,  although  they 
were  for  smaller  scale  earthquakes  [Ref.  2],  Few 
seismologists  now  doubt  that  physical  precursors  to 
earthquakes  do  in  fact  exist.  However,  the  ability  to  be 
sufficiently  consistent  and  uniform  enough  to  permit 
development  of  a  routine  and  reliable  prediction  system 
does  not  exist. 

B.  RESEARCH  PROJECT  BACKGROUND 

This  thesis  advances  the  previous  research  into  the 
evidence  of  radio  frequency  (RF)  emissions  as  a  precursor 
to  earthquake  activity.  Citizen  band  and  amateur  radio 
operators  in  the  Hollister,  California  area  have  reported 
an  "increase  in  the  background  noise  level  preceding 
earthquake  activity"  at  27  MHz  and  at  the  six  and  two 
meter  bands.  The  increase  in  noise  was  on  the  order  of 
10  dB  and  preceeded  a  quake  by  12  to  24  hours.  At  one 
hour  to  fifteen  minutes  prior  to  an  actual  earthquake, 
noise  levels  returned  to  normal  [Ref.  3:p.  11].  The 
appearance  of  a  phenomena  known  as  earthquake  lights 
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[Ref.  4]  and  the  reception  of  18  MHz  radio  noise  [Ref. 

5]  are  further  examples  of  the  potential  link  between  RF 
energy  and  the  occurrence  of  earthquakes.  Further 
presentation  of  measured  precursory  events  and  the 
associated  theories  to  explain  them  are  presented  later 
in  this  thesis. 

The  design  of  a  freestanding  VHF  noise  measurement 
and  recording  system  to  predict  the  intensity  and 
location  of  earthquakes  from  electromagnetic  radiation  at 
an  active  geophysical  location  along  the  San  Andreas  Rift 
Zone  (fault)  was  initiated  in  1984  [Ref.  6]  by  a  joint 
project  of  the  U.S.  Geological  Survey  Office  and  the 
Naval  Postgraduate  School.  The  following  year  the 
proposed  design  of  a  space  based  sensor  to  predict  the 
intensity  and  location  of  earthquakes  from 
electromagnetic  emissions  was  completed  and  the  initial 
installation  of  the  noise  measurement  and  recording 
system  was  made  [Ref.  3].  It  is  at  this  point  that  the 
endeavors  of  this  thesis  begin. 

The  data  acquired  from  the  VHF  noise  measurement 
system  was  observed  for  background  noise  level 
variations.  Se  i  smo  1  o  g  ic  a  .1  data  acquired  from  the 
Geological  Survey  Office  was  compared  to  the  data 
obtained  on  the  chart  recorder  paper.  In  order  to 
accurately  correlate  the  data  it  was  necessary  to 
manually  measure  the  magnitude  of  the  VHF  signal  at 
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predetermined  time  intervals  off  of  the  chart  recordings, 
then  enter  these  readings  into  a  computer  for  future 
graphical  display.  This  process  of  digitizing  the  data 
every  two  weeks  for  the  recordings  of  ten  separate 
channels  of  information  stored  on  very  lengthy  rolls  of 
chart  recorder  paper  is  extremely  time  consuming.  An 
example  of  some  of  the  first  data  retrieved  from  the 
chart  paper  is  shown  in  Figure  1.1.  "ach  of  the  plots 
represent  two  weeks  of  data  for  one  specific  channel. 
Approximately  six  man  hours  were  used  to  obtain  these 
graphs . 

A  solution  to  reducing  the  amount  of  time  required  to 
analyze  the  collected  data,  as  well  as  obtaining  a  means 
of  transmitting  this  data  over  a  digital  satellite  link 
was  being  sought  after  by  the  U.S.  Geological  Survey 
Office.  A  monitoring  system  in  the  Parkfield,  California 
area  relying  on  the  phone  system  for  the  transmission  of 
data  had  lost  vital  information  during  an  actual 
earthquake. 

C.  FOCUS  OF  STUDY 

In  order  to  fulfill  the  need  for  an  accurate  and 
reliable  radio  frequency  noise  measurement  and  recording 
system  for  detection  on  the  30  MHz  and  150  MHz  bands 
which  provides  the  means  for  rapid  data  acquisition  and 
analysis,  the  design  of  a  computer  controlled  monitoring 
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system  is  presented.  Initially  the  data  is  stored  on 
magnetic  disk  to  allow  for  a  series  of  software  routines 
to  graphically  display  the  data.  In  addition,  the  need 
for  a  survivable  remote  sensing  capability  is  fulfilled 
by  the  design  of  digital  satellite  link  interface.  The 
data  may  be  sent  from  the  constructed  computer  monitoring 
system  to  a  Data  Collection  Platform  Radio  Set 
(D.C.P.R.S.)  for  transmission  via  the  Geostationary 
Operational  Environmental  Satellite  (G.O.E.S.) 
transponder  to  the  Command  and  Data  Acquisition  Station 
in  Wallops  Station,  Virginia.  The  establishment  of  a 
correlation  between  increased  background  RF  noise  to 
future  earthquake  fault  activity  may  be  accomplished  as  a 
result  of  the  improved  data  analysis  capability. 
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II.  BACKGROUND 


A.  STRUCTURES  OF  EARTHQUAKE  PREDICTION 

A  successful  prediction  must  include  the  correct 
assessment  of  three  elementary  factors,  namely  the  time, 
place  and  size  of  the  predicted  earthquake.  In  order  to 
approach  thi3  goal,  a  prediction  should  develop 
systematically  through  the  stages  shown  in  Figure  2.1 
[Ref.  7:p.  195].  The  size  of  the  various  rectangles 

shown  in  the  diagram  represent  the  relative  degrees  of 
uncertainty  in  the  different  types  of  prediction.  A 
discussion  of  the  various  realms  of  the  prediction 
process  follows. 

l.  St.aLlatlfi.al  Prediction 

The  first  stage  of  earthquake  prediction  called 
statistical  prediction,  is  based  on  the  assumption  that 
earthquakes  occur  in  a  sequence  with  a  statistical 
characteristic  which  does  r.ot  change  with  time.  If  the 
sequence  of  the  earthquakes  has  either  a  predominant 
period  or  a  probable  correlation  to  a  known  external 
factor,  then  the  future  activity  may  be  predicted 
statistically.  However,  the  prediction  will  involve  a 
considerable  degree  of  uncertainty.  This  type  of  broad 
scale  prediction  is  satisfactory  for  preliminary  work 
such  as  in  long-term  planning  of  optimum  monitoring 
systems  and  disaster  prevention  work, [Ref.  7:p.  194] 
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2.  .  Tectonic  Prediction 

This  type  of  prediction  is  primarily  concerned 
with  the  magnitude,  type  and  other  tectonic  parameters  of 
an  earthquake  that  is  likely  to  occur  in  a  given  locality 
[Ref.  8:p.  1082],  The  accumulation  law  of 
seismotectonic  movement  relates  the  amount  of  accumulated 
strain  and  stress  across  a  fault  to  the  potential  for  the 
occurrence  of  an  earthquake.  This  type  of  prediction 
will  provide  a  mean  lapse  time  for  the  seismic  cycle,  but 
precise  time  information  is  not  available  [Ref.  7 :  p . 
1951. 

3.  Physical  Prediction 

Physical  prediction  attempts  to  precisely 
determine  the  three  critical  factors  of  time,  place  and 
size  by  the  recognition  of  meaningful  seismic  precursors 
to  major  earthquakes.  it  is  an  essential  requirement 
that  a  sufficient  understanding  of  the  physical  laws  of 
earthquakes  be  available  and  precise  monitoring  of  the 
area  in  and  around  an  earthquake  be  made  at  the  proper 
time.  Preseisrr.ic  ananolies  have  been  noticed  in  the 
study  of  many  other  geophysical  phenomena  [Ref.  9].  In 
addition  to  the  extraordinary  behavior  of  fish  and 
animals,  there  is  an  extensive  set  of  precursors 
experienced  at  various  times  and  in  various  regions  of 
the  world  of  the  following  types: 
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land  deformation 


-  tilt  and  strain 

-  foreshock 

-  microseismicity 

-  source  mechanism 

-  fault  creep  anomaly 

-  seismic  wave  velocity 

-  geomagnetic  field 

-  telluric  (earth)  currents 

-  electromagnetic  emissions 

-  resistivity  of  the  earth 

-  radon  content 

-  underground  water 

-  oil  flow 

An  in  depth  discussion  of  the  shock  wave  theory 
of  radio  emissions  from  earthquake  fault  lines  is 
presented,  as  this  is  the  focal  point  of  the  thesis.  The 
design  and  construction  of  the  computer  controlled  VHF 
radio  noise  measurement  and  recording  system  is  based  on 
the  concept  of  using  electromagnetic  emissions  as  a 
precursor  to  future  earthquake  activity. 

4.  Quick  Alarm 

In  strictest  terms  this  last  stage  is  not  really 
prediction,  but  actually  a  method  of  practical  disaster 
reduction.  This  stage  relies  upon  a  fully-automated 
monitoring  system  which  is  linked  to  critical  public 
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and/or  industrial  facilities  so  that  they  may  be 
controlled  promptly  in  an  emergency.  .  If  the  system 
detects  and  recognizes  an  imminent  seismic  event,  the 
critical  facilities  such  as  the  utilities  of  gas  and 
electricity  could  be  secured  before  the  ma:r.  seismic 
disturbance  strikes. [ Ref .  7:p.  196] 

Practical  earthquake  prediction  may  eventually 
be  accomplished  by  analyzing  the  accumulated  predictions 
of  each  of  the  four  stages.  The  successful  use  of  these 
predictions  in  an  orchestrated  effort  to  accurately  and 
reliably  predict  earthquakes  is  in  the  foreseeable 
future. 

B.  RADIO  FREQUENCY  EMISSIONS 

This  thesis  investigates  the  potential  of  using 
electromagnetic  emissions  in  the  radio  frequency  range  as 
a  precursor  to  earthquake  activity.  Radio  frequency 
emissions  from  active  earthquake  fault  lines  have  been 
reported  by  radio  amateurs  in  the  HF  and  VHF  band  of 
frequencies.  The  emissions  are  theorized  to  be  caused  by 
the  increasing  mechanical  pressure  at  the  fault  causing 
sheer  stresses  that  convert  the  insulating  rock  into  a 
conductor  (or  semi-conductor)  by  reducing  the  energy  gap 
between  the  valence  and  conduction  bands.  During  the 
stress  increase,  microfractures  cause  shock  waves  that 
accelerate  the  electrons  in  the  conducting  bands  to 
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radiate  both  acoustic  and  electromagnetic  energy.  If  the 
mechanial  forces  increase  to  the  yield  point,  the  rock 
fails  in  brittle  fracture  and  causes  a  sudden  decrease  in 
the  shear  forces  resulting  in  electron  accelerations 
that  produce  either  large  increases  in  radio  noise,  or  in 
the  phenomena  known  as  clear  weather  1 ightning. [ Ref .  10: 

p.  1] 

The  initial  noise  measurement  and  recording  system 
using  chart  recorders  to  amass  data,  was  placed  into 
operation  in  November  of  1985  [Ref.  3:p.  67].  The  data 
on  all  ten  channels  was  being  collected  and  the  chart 
recorder  paper  being  changed  every  two  weeks  as  designed. 
On  January  26,  1986  a  mild  earthquake  occurred  in  the 
general  vicinity  of  the  Hollister  site.  Less  than  1  km 
away  from  the  detection  system  a  local  winery  had  a  large 
cask  of  wine  fall  from  its  stand  as  a  result.  Analysis 
of  the  data  obtained  from  the  chart  recordings  for  the 
time  frame  surrounding  the  incident  appears  to  support 
the  observations  of  the  amateur  and  citizen  band  radio 
operators.  An  increase  in  the  general  background  noise 
floor  prior  to  the  quake  occurred,  as  well  as  a  sharp 
decrease  in  the  level  just  prior  to  the  quake.  Figure 
2.2  illustrates  the  radio  frequency  noise  levels  lor  the 
January  25  to  27 »  1986  period  on  both  the  38  MHz  and  150 
MHz  frequencies. 
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Figure  2.2  Collected  Data  for  25  to  27  January  1986 
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III.  THE  SATELLITE  .LINK 


A.  GENERAL  BACKGROUND 

To  fulfill  the  requirement  that  the  currently 
designed  computer  controlled  VHF  noise  measurement  and 
recording  system  be  capable  of  having  the  collected  data 
transmitted  over  a  satellite  link,  the  U.S.  Geological 
Survey  Office  proposed  the  use  of  a  commercial 
transmitter  which  utilizes  the  Geostationary  Operational 
Environmental  Satellite  Data  Collection  System. 
Commercial  UHF  satellite  traTsmitters  are  currently  being 
used  to  provide  a  remote  sensing  capability  for  ^-j.ected 
environmental  sensors  and  are  compatible  with  this 
system,  [Ref.  11]  and  [Ref.  12]  provide  the 

background  information  for  interfacing  the  computer 
controlled  VHF  r.oise  measurement  and  recording  system  to 
the  established  data  collection  satellite  link.  The 
system  becomes  a  remote  sensing  device,  providing  ready 
access  to  the  data  without  traveling  to  the  site  to 
retrieve  the  data  stored  on  the  magnetic  disk.  Specifics 
of  the  actual  hardware  interfacing  of  the  designed  system 
to  the  satellite  link  are  provided  later  in  the  thesis. 
To  provide  a  workable  understanding  of  the  system  so 
future  users  are  familiar  with  its  overall  operation,  a 
basic  overview  of  the  satellite  link  is  presented  [Ref. 
ll:pp.  2-17]  and  [Ref.  12:pp.  2-9]. 
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There  are  currently  two  functioning  Geostationary 
Onerational  Environmental  Satellites  (G.O.E.S.)  that  are 
part  of  a  system  of  earth  and  space  environmental  sensors 
which  provide  an  almost  continuous  observational 
information  flow  to  the  ground  based  user  stations.  The 
GOES  system  is  operated  and  controlled  by  the  National 
Environmental  Satellite  Data  and  Information  Service 
(NESDIS).  The  organization  was  formerly  called  the 
National  Earth  Satellite  Service  (MESS).  NESDIS  is  a 
suborganization  of  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  which  is  controlled  by  the  U.S. 
Department  of  Commerce.  The  GOES  system  itslf  was  a 
result  of  the  combined  efforts  of  NESDIS  and  the  National 
Aeronautics  and  Space  Administration  (NASA)  based  on  the. 
experiments  with  the  NASA  advanced  technology  satellites. 

The  GOES  system  consists  of  several  subsystems,  one 
of  which  is  the  Data  Collection  System  (DCS).  The  Data 
Collection  System  uses  the  GOES  spacecraft  for  the 
relay  of  data  from  remotely  located  sites  on,  or  nearly 
on  the  surface  of  the  earth.  Of  course  these  remote 
sites  must  be  in  radio  view  of  the  satellite.  There  are 
many  other  uses  for  the  GOES  satellite,  but  the 
discussion  of  these  is  beyond  the  scope  of  this  thesis. 

There  are  two  GOES  satellites  located  in 
geosynchronous  orbits  over  the  equator  at  75  degrees  West 
and  135  degrees  West  longitudes.  The  radio  view  of  the 
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two  spacecraft  is  from  approximately  0  degrees  westward 
to  165  degrees  East  longitude.  There  is  a  limitation  of 
the  coverage  to  approximately  77  degrees  North  ana  South 
latitudes.  See  Figure  3.1.  In  case  of  a  failure  cf 
either  the  "East  "  or  "West"  satellites  just  mentioned, 
there  is  a  third  geostationary  satellite  in  a 
geosynchronous  orbit  halfway  between  them  at  107  degrees 
West  longitude.  This  craft  is  also  used  to  prevent  data 
loss  as  a  result  of  either  satellite  being  temporarily 
out  of  commission  as  a  result  of  an  eclipse  period. 

B.  THE  GOES  DATA  COLLECTION  SYSTEM 

The  GOES  Data  Collection  System  (DCS)  is  a  satellite 
based  system  for  the  collection  of  a  large  variety  of 
environmental  data  obtained  from  almost  any  point  in  the 
Western  Hemisphere.  It  is  a  data  relay  network  composed 
of  literally  thousands  of  individual  data  gathering 
devices  known  as  Data  Collection  Platforms  that  each 
transmit  the  data  they  have  acquired  to  one  of  the  GOES 
satellites.  The  signal  is  then  relayed  from  the 
satellite  to  the  Command  and  Data  Acquisition  (CDA) 
Ground  Station  located  in  Wallops  Station,  Virginia. 

The  transponders  carried  onboard  the  GOES  spacecrafts 
relay  the  UHF  transmissions  from  the  Data  Collection 
Platforms  (ranging  rom  401.7  MHz  up  to  402.1  MHz)  via  S- 
band  transmission  to  the  ground  receiving  station. 
Conversely,  S-band  transmissions  from  the 
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NESDIS  ground  station  can  be  relayed  through  the 
spacecraft  transponder  and  converted  into  a  UHF 
transmission  for  dissemination  to  properly  equipped 
receivers  in  radio  view  of  the  satellite. 

The  satellite  transponders  are  capable  of  supporting 
up  to  233  reply  channels.  There  are  200  channels 
reserved  for  use  by  regional  and  domestic  channels 
ranging  from  401.7  MHz  up  to  402.0  MHz,  each  with  a  1.5 
KHz  channel  separation.  The  other  33  channels  are  for 
international  use  in  the  402.0  MHz  to  402.1  MHz  range. 
Each  has  a  3.0  KHz  channel  separation.  The  ground  system 
is  limited  to  supporting  80  channels  and  5000  Data 
Collection  Platforms. 

There  are  four  functional  subsystems  of  the  GOES  Data 
Collection  System.  They  are  as  follows: 

-  Deployed  Data  Collection  Platform  Radio  Sets 

-  East  and  West  Spacecraft  Transponders 

-  Command  and  Data  Acquisition  Station 

-  Central  Data  and  Dissemination  Facility 

A  brief  discussion  of  each  of  these  is  necessary  for  a 
basic  understanding  of  how  the  GOES  Data  Collection 
System  operates  and  subsequently  how  the  designed  VHF 
noise  measurement  and  recording  system  can  interface  with 
the  satellite  link  to  achieve  the  remote  sensing 
capability.  Figure  3.2  shows  a  block  diagram  of  the 
G.O.E.S.  Data  Collection  System. 
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1.  Data  Collection  Platform  Radio  Sets 


There  are  five  different  types  of  Data 
Collection  Platform  Radio  Sets  (DCPRS)  which  permit 
various  modes  of  system  operation.  They  are  as  follows: 

-  Self-Tined 

-  Self-Timed  and  Random  Reporting 

-  Random  Reporting 

-  Interrogated 

-  International  Self-Timed 

A  description  of  each  type  follows, 
a.  Self-Timed 

Self-Timed  DCPRS  contain  only  a  transmitter 
and  a  timing  source  which  is  preprogrammed  by  the  user  to 
report  during  a  specific  hour  at  a  specific  minute  and  at 
a  specific  rate  over  a  ZH  hour  period.  The  preprogrammed 
reporting  time  and  rate  of  the  DCPRS  is  not  at  all 
related  to  the  cycle  of  the  sensor  attached  to  the 
platform.  The  sensor  cycle  refers  to  the  time  interval 
or  frequency  of  when  the  environmental  sensor  attached 
directly  to  the  DCPRS  reads  its  specific  data.  The  data 
is  then  placed  into  a  memory  buffer  in  the  DCPRS  for 
transmission  at  the  preprogrammed  time.  For  example,  a 
DCPRS  may  have  data  transferred  to  it  from  an 
environmental  sensor  every  five  minutes,  but  this  data 
would  be  stored  temporarily  and  then  transmitted  through 
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the  GOES  Data  Collection  System  only  once  every  three 
hours . 

b.  Self-Timed  and  Random  Reporting 

Self-Timed  and  Random  Reporting  DCPRS  have 

the  same  characteristics  as  the  previously  mentioned 
type,  however,  it  also  has  the  ability  to  transmit  over  a 
secondary  channel  when  environmental  conditions  require 
more  frequent  reporting  than  offered  under  self-timed 
operation.  In  this  DCPRS  the  normal  self-timed  mode  can 
be  used  under  most  circumstances,  but  if  the  need  for 
more  immediate  data  is  required,  the  DCPRS  car.  be 
automatically  placed  into  the  random  reporting  mode. 

c.  Random  Reporting 

The  Random  Reporting  DCPRS  contains  a 
transmitter  which  will  broadcast  at  a  random  time.  This 
is  controlled  by  environmental  conditions  being  monitored 
by  the  sensor  having  reached  a  preassigned  threshold. 
However,  the  Random  Reporting  DCPRS  are  also  preassigned 
to  report  at  least  three  times  a  day  in  order  to  assure 
that  the  DCPRS  and  the  sensor  are  functioning  properly. 

d.  Interrogated 

Unlike  the  other  types  of  DCPRS,  the 
Interrogated  DCPRS  contain  both  a  transmitter  and  a 
receiver.  The  onboard  receiver  is  preset  to  either  the 
East  or  West,  spacecraft  DCPRS  interrogation  link 
frequency  which  enables  it  to  detect  its  own  unique 
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identification  or  address  when  transmitted  from  the 
Command  and  Data  Acquisition  (CDA)  Station  at  Wallops 
Station,  Virginia.  If  a  DCPRS  recognizes  its  own  special 
identification  or  callsign,  it  is  programmed  to  transmit 
all  of  the  data  it  has  stored  since  the  previous 
interrogation.  The  scheduling  of  the  Interrogated  DCPHS 
is  very  similar  to  the  Self-Timed  DCPRS,  except  that  the 
schedule  is  stored  at  the  Command  and  Data  Acquisition 
Station  instead  of  within  the  DCPRS  themselves.  The 
interrogation  of  the  DCPRS  may  occur  at  intervals  as 
small  as  5  minutes  or  as  large  as  once  every  24  hours. 
In  addition,  some  of  the  interrogated  DCPRS  have  a 
secondary  alert  channel  which  will  send  an  alert  signal 
to  the  CDA  telling  it  to  interrogate  because  a 
preassigned  threshold  has  been  reached. 

e.  International  Self-Timed 

The  International  Self-Timed  is  identical  to 
the  domestic  Self-Timed  except  for  the  range  over  which 
the  UHF  transmitters  may  transmit. 

2.  East  and  West  Spacecraft  Transponder? 

The  next  subsystem  of  the  GOES  Data  Collection 
System  to  be  discussed  is  the  set  of  transponders  aboard 
the  East  and  West  spacecraft.  The  GOES  "East"  spacecraft 
is  located  at  75  degrees  West  and  the  GOES  "West" 
spacecraft  at  135  degrees  West.  Each  of  the  craft  are 
capable  of  handling  the  233  channels  allocated  for  the 
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system  (200  domestic  and  33  international).  However, 
extra  channel  bandwidth  is  provided  for  all  of  the 
domestic  channels  by  assigning  the  odd  numbered  channels 
to  the  East  spacecraft  transponder  and  the  even  channels 
to  the  West  spacecraft  transponder.  As  a  result,  a  3.0 
KHz  guard  bandwidth  is  achieved  between  adjacent 
channels.  A  further  organization  of  the  channels  is  made 
by  the  assignment  of  certain  blocks  of  channels  on  each 
spacecraft  to  specific  types  of  DCPRS.  See  Figure  3-3 
for  a  comparison  of  the  odd/even  channel  assignments. 
The  international  channels  are  assigned  to  both  the  East 
and  the  West  spacecraft  at  all  times. 

The  Command  and  Data  Acquisition  Station  sends 
out  command  and  interrogate  signals  to  the  appropriate 
spacecraft  transponder  at  an  S-band  frequency  that  is 
transformed  into  a  UHF  signal.  Each  spacecraft 
transponder  sends  out  a  different  UHF  frequency  to  the 
Data  Collection  Platform  Radio  Sets.  The  assigned 
frequency  for  the  East  is  468.8375  MHz  and  the  assigned 
West  frequency  is  468.825  MHz.  The  response  signals  sent 
up  from  the  Interrogated  DCPRS  and  the  self-transmitted 
signals  from  the  Self-Timed  and  the  Random  Reporting 
DCPRS  are  sent  up  at  UHF  frequencies.  The  signals  are 
then  translated  by  the  satellite  transponders  to  S-band 
and  then  retransmitted  to  the  Command  and  Data 
Acquisition  Station. 


Figure  3«3  East  and  West  Spacecraft  Channel  Numbering 
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3 .  Command  and  Data  Acquisition  Station  (CPA) 


The  messages  are  received  from  either  one  of  the 
spacecraft  via  their  S-band  transmissions  into  the  CDA. 
The  received  Data  Collection  Platform  Racio  Set  messages 
are  then  routed  to  the  Data  Acquisition  and  Monitoring 
Subsystem  (DAMS)  units  which  will  demodulate  the  data  and 
perform  signal  quality  measurements.  A  typical  DCP 
message  containing  the  DAMS  quality  information  is  shown 
in  Figure  3.4.  The  data  information  is  then  transferred 
from  the  CDA  to  the  Central  Data  and  Dissemination 
Facility  over  high  speed  dedicated  circuits  for  further 
dissemination.  As  mentioned  previously,  the  CDA  computer 
contains  the  schedule  files  for  the  Self- Timed  and 
Interrogated  types  of  DCPRS. 

4.  Central  Data  and  D.Lsfi.emiiiaUpn  Facility 

The  Central  Data  and  Dissemination  Facility 
(CDDF)  is  located  in  the  World  Weather  Building  at  Camp 
Springs,  Maryland.  There  are  two  identicial  computers 
there.  One  of  them  is  on-line  at  all  times,  while  the 
other  is  used  for  statistical  analysis.  The  incoming 
messages  from  the  Command  and  Data  Acquisition  Station 
are  then  routed  to  the  users  of  the  DCPRS  via  landline 
along  with  any  pertinent  error  information  detected 
during  the  processing  of  the  message. 
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C.  DIGITAL  SATELLITE  LINK  INTERFACE 

The  designed  computer  controlled  VHF  noise 
measurement  and  recording  system  is  capable  of  sending 
binary  data  out  of  an  RS232  port  in  a  serial  format 
with  a  number  of  software  programmable  variations  of  the 
data  bit  stream.  A  detailed  analysis  of  how  the  system 
can  forward  the  acquired  data  to  the  buffer  in  a  Data 
Collection  Platform  Radio  Set  will  be  presented  later  in 
this  thesis. 
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IV.  DESCRIPTION  OF  RESEARCH 


A.  DESIGN  OVERVIEW 

The  design,  construction  and  implementation  of  the 
computer  controlled  radio  frequency  (RF)  noise 
measurement  and  recording  system  was  accomplished  by 
advancing  the  previous  research  into  precursory 
electromagnetic  activity  associated  with  earthquakes 
[Ref.  3]  and  [Ref.  6],  The  location  of  the  installed 
monitoring  system  is  three  miles  east  of  San  Juan 
Batista,  California  and  seven  miles  south  of  Hollister. 
This  site  is  located  within  the  Hollister  Hills  State 
Vehicular  Recreation  Area  (upper  ranch),  directly  south 
of  the  MX  track.  The  site  was  chosen  in  previous 
research  projects  because  it  is  historically  an  active 
location.  In  addition,  the  location  is  remote  enough  to 
be  considered  a  rural  area  with  respect  to  man-made  radio 
noise.  Figure  4.1  refers.  The  exact  location  is  one 
thousand  feet  from  the  San  Andreas  fault  at  121  degrees 
and  23.5  minutes  West  longitude  and  36  degrees  45.5 
minutes  North  latitude.  Figures  4.2  and  4.3  provide  a 
detailed  means  of  locating  the  site. 

The  completion  of  the  thesis  involved  the 
coordination  of  efforts  in  both  hardware  and  software 
design  to  accomplish  the  established  goals  of  obtaining 
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a  data  recording  and  analysis  system,  as  well  as  the 
digital  satellite  link.  Figure  4.4  provides  a  block 
diagram  of  the  hardware  system  and  a  flow  chart  of  the 
software  used  on  the  Maval  Postgraduate  School  IBM  3C 3 3 
computer  to  provide  graphical  display  of  the  data  stored 
on  5.25  inch  magnetic  disks.  The  two  frequencies  being 
monitored  (38.45  MHz  and  150.75  MHz)  have  a  bandwidth  of 
50  KHz.  They  have  been  established  as  receive  only 
channels  to  allow  the  Naval  Postgraduate  School  to 
investigate  the  short  and  long  term  variations  of 
backround  noise  (noise  floor).  The  Army  Frequency 
Coordinator  for  the  Western  United  States  will  continue 
these  assignments  on  a  yearly  basis. 

The  antenna  structure  consists  of  two  twenty  foot 
towers  and  ten  antennas.  Each  tower  has  five  antennas 
mounted  on  it.  Eight  yagi  antenna  receive  the 
horizontally  and  vertically  polarized  signals  of  the  two 
frequency  ranges.  Omni  antenna  for  each  of  the 
frequencies  are  vertically  mounted.  See  Figure  4.5  for  a 
list  of  the  specific  channel  assignments.  The  signals 
received  by  each  of  the  antenna  are  individually  cabled 
to  ten  modular  amplitude  modulation  (AM)  receivers. 

Once  the  receiver  has  processed  a  signal,  a  voltage 
level  corresponding  to  the  strength  of  this  incoming 
signal  is  obtained  from  the  "S"  meter  output.  This 
output  voltage  level  is  sent  to  an  interface  amplifier 


Figure  4.4  Design  Overview  Diagram 
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CHANNEL  ASSIGNMENTS 

FREQUENCY 


CHANNEL 

(MHZ) 

POLARIZATION 

DIRECTION 

1 

150.75 

VERTICAL 

NORTHWEST 

* 

\ 

i 

150.75 

VERTICAL 

SOUTHEAST 

!  , 

-  -j* 

l 

150.75 

OMNI 

4 

150.75 

VERTICAL 

NORTHWEST 

5 

150.75 

HORIZONTAL 

SOUTHEAST 

6 

33.  45 

VERTICAL 

SOUTHEAST 

7 

38.45 

HORIZONTAL 

SOUTHEAST 

3 

38.  45 

VERTICAL 

NORTHWEST 

9 

38.  45 

HORIZONTAL 

NORTHWEST 

10 

33.45 

OMNI 

Figure  4.5 

Channel  Assignments 
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where  it  is  amplified  to  a  range  of  zero  to  five  volts, 
corresponding  to  the  minimum  and  maximum  received  signals 
of  interest,  respectively.  The  zero  to  five  volt  range 
signal  is  then  quantized  by  an  8-bit  analog  to  digital 
converter.  The  data  is  stored  in  an  ASCII  format  on 
magnetic  disk.  The  data  is  also  sent  out  RS232  ports  to 
a  buffer  on  a  Self-Timed  Data  Collection  Platform  Radio 
Set  (D.C.P.R.S.)  for  later  transmission  over  a  satellite 
link,  as  well  as  to  a  data  monitor  for  visual  display  of 
the  sampled  signal  level  data.  Control  of  this  process 
is  obtained  from  a  program  stored  in  an  erasable 
programmable  read  only  memory  (EPROM)  chip  on  the  single 
board  computer.  Power  for  the  system  is  obtained  from 
either  an  AC  or  DC  source. 

Data  analysis  is  accomplished  by  plotting  the  data 
stored  on  magnetic  disk  by  using  the  IBM  3033  network  at 
the  Naval  Postgraduate  School.  The  data  is  stored  in  a 
very  compact  manner  onto  a  formatted  disk.  However,  the 
disk  must  be  retrieved  once  a  month  and  a  new  one 
inserted.  At  this  point  the  series  of  software  programs 
illustrated  in  Figure  4.4  are  performed  to  process  the 
data.  The  raw  aata  is  first  rearranged  into  a  PCDOS  file 
format.  The  PCDOS  file  is  then  transferred  to  the  main 
frame  using  an  established  microcomputer  link.  The  data 
is  arranged  into  a  sequential  format  by  the  next  program. 
In  the  final  stage  of  the  processing  the  graphical  plot 
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is  obtained.  The  complete  process  for  obtaining;  a 
graphical  representation  of  thirty  days  of  signal  level 
information  for  ten  separate  channels  has  been  reduced  to 
hours,  as  opposed  to  days  or  even  weeks  of  effort 
required  prior  to  implementation  of  this  system.  The 
collected  data  is  easily  and  rapidly  correlated  to  any 
additional  earthquake  related  data  provided  oy  the 
Geological  Survey  Office  in  Menlo  Park,  California. 

B.  SYSTEM  HARDWARE  DESIGN  AND  OPERATION 

The  block  diagram  of  the  hardware  components  of  the 
computer  controlled  RF  noise  measurement  and  recording 
system  is  presented  again  in  greater  detail  in  Figure 
4.6.  A  close  analysis  of  each  of  the  individual 
subsystems  is  presented, 
l .  Efca&ltt&na 

Each  one  of  the  Hamtronics  VHF  converter  modules 
are  designed  to  amplify  and  convert  a  specific  frequency 
to  the  HF  range.  Schematic  diagrams  of  the  38.45  MHz  and 
the  150.75  MHz  converters  are  shown  in  Figure  4.7  and 
Figure  4.8,  respectively.  Each  one  of  the  converters 
have  10.7  MHz  intermediate  frequency  (IF)  transformers 
for  the  mixer  output  circuits.  Tha  mixer  output  circuits 
themselves  consist  of  a  slug  tuned  coil,  capacitive 
voltage  dividers  and  an  RCA  jack  to  provide  a  50  ohm 
output. 
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The  10.7  MHz  IF  signal  is  sent  from  tne 
converter  module  to  the  IF/Audio  module  shown  in  Figure 
4.9.  The  IF/Audio  assembly  is  a  highly  sensitive  and 
selective  IF  ampliier,  AH  detector,  audio  amplifier  and 
squelch  system  which  is  used  in  conjunction  with  the 
converter  to  create  a  VHF  communications  receiver.  The 
HF  signal  is  then  mixed  with  the  10.245  MHz  oscillator 
frequency  and  then  down-shifted  to  455  KHz.  This  signal 
is  then  amplified  and  detected.  The  MS"  meter  output  to 
the  interface  amplifier  is  obtained  from  an  integral  log 
detector. 

The  selectivity  of  this  unit  (plus  or  minus  10 
KHz  at  -6dB  and  plus  or  minus  25  KHz  at  -60  dB)  will 
permit  its  use  as  a  noise  floor  measurement  system.  An 
optimally  coupled  three  stage  IF  filter  provides  the 
selectivity  against  adjacent  channel  interference.  The 
unit  has  automatic  gain  control  (AGO  for  both  IF  and 
RF,  with  the  RF  AGO  being  delayed  for  best  response. 
Each  receiver  is  installed  in  its  own  individual  aluminum 
casing  to  prevent  the  detrimental  effects  of 
electromagnetic  inteference  (EMI). 

2.  ;nterfacs,  MplLf.ier^ 

The  "S"  meter  output  from  the  integral  log 
detector  in  the  receiver  is  sent  directly  to  the 
interface  amplifier  assembly  shown  in  Figure  4.10. 
Previously  the  "S"  meter  output  and  a  time  pulse  from  a 
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Figure  4.9  IF/Audio  Assembly 


timer  assembly  were  integrated  into  one  signal  [Ref.  3: 
p.  58].  Observation  of  the  data  found  the  design  to  be 
sensitive  to  temperature  and  subject  to  nonlinear 
variations  in  output  levels.  In  the  current  system  a 
great  deal  of  simplification  has  resulted  by  utilizing 
the  computer  as  the  timing  and  control  device.  The 
interface  amplifier  uses  LM747  series  general  purpose 
dual  operational  amplifiers.  As  a  result,  no  frequency 
compensation  is  required  and  there  is  an  inherent 
balanced  offset  null.  The  "S"  meter  output  of  the 
receiver  is  amplified  by  the  first  stage,  then  inverted 
to  a  positive  voltage  using  a  unity  gain  in  the  second 
stage. 

Previous  data  monitoring  indicates  that  the 
received  input  signals  fell  in  the  range  of  zerc  to  ten 
microvolts.  Amplification  of  these  signals  to  the 
millivolt  range  is  available  at  the  "S"  meter  output  of 
the  receivers.  The  signal  is  then  amplified  to  the 
required  range  of  zero  to  five  volts  by  the  interface 
amplifier  assembly,  correspond iong  to  the  required  input 
range  of  the  analog  to  digital  converter.  Figure  4.11 
provides  a  graphical  plot  of  the  input  RF  signal  to  the 
receiver  versus  the  corresponding  analog  to  digital 
converter  input. 
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RCVR  INPUT  VS  A/D  CONV  INPUT 


(SXIOAOOH)  lAdN'I  H3A13D3H 


Figure  4.11  Graph  of  Receiver  Input  vs 
A/D  Converter  Input 
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«.0  S.0  4.0  0.0 

A/D  CONVERTER  INPUT  (VOLTS) 


3.  Analog  To  Digital  Converter 


Quantizing  of  the  zero  to  five  volt  signal  level 
from  the  interface  amplifier  assembly  is  accomplishec  by 
an  8-bit  analog  to  digital  converter.  The  16  channel 
8-fcit  A/D  converter  with  the  associated  address  decoder 
and  multiplexing  circuitry  are  located  on  the  Vesta 
OEM188  multifunction  board  shown  in  Figure  4.12.  The 
board  has  been  jumpered  to  allow  single-ended  analog 
inputs.  Port  5  (P5)  on  the  multifunction  board  will 
allow  for  up  to  16  channels  of  single-ended  analog  input. 
(P5)  is  set  up  as  shown  in  Table  I. 

The  multifunction  card  has  a  bus  interface  to 
allow  for  direct  connection  to  the  microprocessor  board. 
The  bus  interface  has  a  wait-state  generator  to  induce  a 
1.2  microsecond  hold  state  to  stabilize  the  bus  during 
input  and  output  procedures.  The  wait-state  generator 
will  slow  down  the  microprocessor  to  allow  for  the 
operation  of  slow  input/output  devices  such  as  the  A/D 
converter.  The  bus  is  configured  as  a  60  conductor,  dual 
row  header  that  is  designated  Port  1  CPI)  on  the 
multifunction  card.  The  header  is  connected  via  ribbon 
cable  to  the  expansion  bus  connector  on  the  single  board 
computer  itself. 

4.  Computer 

The  Vesta  OEM188  single  board  computer  using  an 
Intel  8Cuu  „;icroprocessor  with  a  clockspeed  of  8  MHz 
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TABLE  I 


MULTIFUNCTION  CARD  PORT  5  (P5)  PIN  ASSIGNMENT 
(A/D  CONVERTER  INPUTS) 

HBER  FUNCTION  ASSIGNMENT 


22 

Analog 

(  +  ) 

0 

Channel  1 

input 

24 

Analog 

(  +  ) 

1 

Channel  2 

Input 

26 

Analog 

<♦> 

2 

Channel  3 

Input 

16 

Analog 

(  +  ) 

3 

Channel  4 

Input 

23 

Analog 

U> 

4 

Channel  5 

Input 

18 

Analog 

(  +  ) 

5 

Channel  6 

Input 

25 

Analog 

(  +  ) 

6 

Channel  7 

Input 

20 

Analog 

<♦> 

7 

Channel  8 

Input 

9 

Analog 

(  +  ) 

8 

Channel  9 

Input 

11 

Analog 

<♦) 

9 

Channel  10 

Input 

13 

Analog 

(♦) 

10 

Unassigned 

19 

Analog 

<+) 

11 

Unassigned 

15 

Analog 

<♦) 

12 

Unassigned 

21 

Analog 

<+) 

13 

Unassigned 

17 

Analog 

<♦> 

14 

Unassigned 

7 

Analog 

U) 

15 

Unassigned 

8 

Ground 

Ground 

10 

+  5  DC 

+  5  DC 
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provides  a  maximum  level  of  control  with  minimal  size 
requirements.  The  control  program  for  the  microprocessor 
has  been  stored  on  an  EPROM,  allowing  the  disk  drives 
attached  to  the  computer  to  be  used  solely  for  data 
storage.  Both  the  operating  system  basic  input/output 
system  (BIOS)  and  higher  level  language  compiler  (BASIC) 
are  also  stored  in  ROM  directly  on  the  board.  The 
computer  has  a  battery-backed  real  time  clock  (National 
Semiconductor  MM 5 827 4)  mounted  on  it,  as  well  as  a  Y/D1770 
disk  drive  controller  which  will  support  up  to  four  5.25. 
inch  disk  drives.  In  addition,  an  EPROM  programmer  is 
located  in  socket  U 3 1  on  the  board  to  allow  for  direct 
storage  of  programs  written  on  the  computer  in  BASIC.  A 
schematic  of  the  computer  is  provided. in  Figure  4.13. 

The  terminal  connector  located  at  Port  3  (P3)  is 
designed  to  drive  a  high  speed  full  duplex  CRT  terminal. 
This  is  the  primary  means  of  operator  communications  with 
the  computer.  The  three  pins  on  the  P3  connection  and 
their  corresponding  functions  are  described  in  Table  II. 
Communications  take  place  over  an  RS-232  line  at  a  rate 
selected  by  the  configuration  switch  located  on  the 
single  board  computer  itself.  The  function  of  this 
switch  is  twofold.  The  switch  provides  the  system  with 
information  concerning  the  extent  of  low  memory.  The 
switch  also  determines  the  baud  rate  at  which  the  CRT 
serial  channel  will  communicate.  Table  III  provides  the 
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Figure  4.13  Schematic  Diagram  of  OEM188  Computer 
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TABLE  II 


OEM188  PORT  3  (P3)  PIN  ASSIGNMENTS 
(TERMINAL  CONNECTOR) 


TXD,  transmitted  data  from  computer 

GND,  signal  common 

RXD,  received  data  from  terminal 


TABLE  III 

OEM  188  CONFIGURATION  SWITCH  SETTINGS 

-&EIIIHfiS  F-U-NCIIP.M 


0  X  X  X 


1  X  X  X 


<  64K  Onboard  low  frequency 
>  128K  Onboard  low  memory 


X  0  0  0 


150  BAUD 


X  0  0  1 


300  BAUD 


X  0  1  0 


600  BAUD 


X  0  1  1 


1200  BAUD 


X  1  0  0 


2400  BAUD 


X  1  0  1 


4800  BAUD 


X  1  1  0 


9600  BAUD 


X  1  1  1 


19200  BAUD 
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correct  settings  for  this  switch.  Switch  positions  are 
shown  from  left  to  right  corresponding  to  pin  numbers  1 
to  4  (  0  =  closed,  1  =  open).  The  current  system 
incorporates  only  1 6k  of  onboard  low  memory.  The  present 
configuration  switch  setting  is  0-1-1-0  to  accommodate 
terminal  interfacing  at  a  9600  baud  rate. 

The  power  connector  located  at  port  1  (PI) 
supplies  the  system  with  +5  VDC  and  +12  VDC.  Pin  1  of  PI 
is  connected  to  +5  VDC  and  pin  3  connected  to  ground. 
The  +12  VDC  connection  is  not  required.  As  stated 
previously,  up  to  four  5.25  inch  floppy  disk  drives  may 
be  used.  In  the  present  arrangement,  Port  7  ( P7 )  has 
been  connected  to  two  disk  drives  via  34  conductor  ribbon 
cable  in  a  daisy  chain.  Power  is  supplied  to  the  disk 
drives  separately  from  the  system's  power  supply. 

Auxiliary  RS-232  serial  port  4  ( P 4 )  on  the 
computer  board  is  used  to  transmit  the  collected  data  to 
the  Data  Collection  Platform  Radio  Set  buffer  for 
transmission  over  the  satellite  link  at  a  later  time. 
The  information  is  sent  out  this  port  in  a  software 
programmable  format.  Table  IV  describes  the  pin 
assignments  and  their  functions. 

5.  Disk  Drives 

Two  Teac  model  FD-55B  half-height  double-sided 
double-density  disk  drives  are  installed  in  the  system. 
The  data  disk  is  installed  in  the  "A"  drive  to  receive 
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AUXILIARY 

TABLE  IV 

RS-232  SERIAL  PORT  4  (P4)  PIN 

ASSIGNMENTS 

PIH  HUMBER 

FUNCTION 

1 

TXD ,  transmitted  data 

from  computer 

2 

GND,  signal  common 

3 

RXD,  received  data 

4 

RTS,  ready  to  3end 

5 

CTS,  clear  to  send 

_ 
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the  digitized  signal  levels  of  all  ten  channels  at  each 
sampling  interval.  The  "B"  drive  was  used  in  the  initial 
design  of  the  computer  control  program  and  may  be 
alternated  with  the  currert  "A"  drive  for  data  storage. 
Both  drives  have  been  sealed  in  a  dust,  insect  and  vermin 
proof  plexiglass  container  to  prevent  any  possible  damage 
to  moving  parts,  the  writing  heads,  or  the  data  disk 
themselves . 

Each  drive  has  a  power  requirement  of  +  12  VDC  at 
0.9  amp  max  and  +  5  VDC  at  0.6  amp  max.  Typical  power 
consumption  while  in  operation  is  3.9  watts  and  while 
idle  2.6  watts.  The  units  are  extremely  reliable,  with 
an  operating  temperature  range  of  4  degrees  to  46  degrees 
centigrade  and  a  relative  humidity  operating  range  of 
205  to  805  (noncondensing). 

6.  Power  Supply.  Sxalsm 

The  system  has  been  designed  to  operate  using 
either  60  Hz,  115  volt  AC  power,  or  12  volt  DC  power. 
Figure  4.14  shows  the  block  diagram  of  the  power  supply 
system.  When  AC  power  is  available,  the  system  is 
designed  to  supply  both  plus  and  minus  12  volts  DC  and 
plus  five  volts  DC.  The  battery  chargers  will  maintain  a 
constant  charge  on  the  12  volt  DC  batteries  until  AC 
power  is  lost.  At  that  time,  the  batteries 
instantaneously  begin  to  provide  the  system  with  the 
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required  voltages.  The  batteries  will  maintain  the 
system  in  full  operation  for  approximately  48  hours. 
Even  if  the  system  were  forced  to  operate  without  AC 
power  for  a  period  of  time  sufficient  enough  to  drain  the 
power  batteries,  the  small  nickel-cadmium  battery  on  the 
computer  board  will  provide  power  to  the  real  time  clock 
for  up  to  a  week.  Although  data  storage  capability  will 
be  lost  during  this  time,  the  automatic  reset  of  the 
computer  will  reactivate  the  control  program  immediately 
once  AC  power  is  restored. 

7.  Interface  Connections 

Two  25  pin  male  D-type  connectors  are  installed 
on  the  front  of  the  system  enclosure  to  provide  serial 
data  output  to  the  external  peripherals.  Output  Port 
number  1  is  wired  to  the  terminal  connector  Port  3  (P3) 
of  the  computer.  The  pir.  connections  and  their 
corresponding  functions  are  available  in  Table  II.  The 
baud  rate  at  which  the  communication  of  data  is  sent  out 
of  output  port  1  is  set  by  the  configuration  switch  on 
the  computer.  Table  III  provides  the  available  options 
and  correct  settings  for  the  switch.  The  switch  is 
presently  set  to  0-1-1-0  to  accommodate  full  duplex 
terminal  interactions  at  a  9600  baud  rate. 

Data  output  port  number  2  is  connected  to 
auxiliary  RS-232  port  4  (P4)  on  the  computer.  At  the 
time  of  each  sampling  interval  the  signal  level  data  is 
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sent  out  of  this  port  in  a  serial  format  to  the  Data 
Collection  Platform  Radio  Set  buffer.  The  data  is  stored 
in  the  buffer  until  it  is  transmitted  over  the  402  MHz 
satellite  channel  at  the  predetermined  interval.  Table 
IV  describes  the  pin  assignments  and  their  functions. 
The  software  design  section  describes  in  detail  the 
options  available  ir.  formatting  the  data  output.  The 
data  is  presently  being  sent  out  at  a  rate  of  110  baud 
from  output  Port  2  to  the  D.C.P.R.S  buffer  as  timing  is 
not  at  all  critical.  The  baud  nte  .an  be  increased  to 
as  high  as  9600  baud  if  necessary. 

C.  SYSTEM  SOFTWARE  DESIGN  AND  OPERATION 

The  flow  chart  of  the  software  programs  used  to 
process  the  recorded  data  of  the  RF  noise  measurement 
system  is  presented  again  in  Figure  4.15.  A  description 
of  each  stage  of  the  data  processing  is  presented. 

1 •  Digitize.com 

The  program  Digitize.com  is  stored  in  the  model 
2764  EPROM  located  on  the  Vesta  OEM 188  single  board 
computer.  Control  of  the  entire  system  is  maintained  by 
this  program.  Appendix  A  contains  a  copy  of  Digitize.com 
written  in  the  higher  level  language  BASIC.  The  computer 
samples  the  ten  analog  input  signal  levels  sent  from  the 
interface  amplifier  assemblies  with  an  8-bit  A/D 
converter  at  a  specific  time  interval.  The  digitized 
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Figure  4.15  Data  Processing  Flow  Chart 
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signal  level  data  is  then  stored  on  magnetic  disk.  In 
addition,  the  data  is  sent  to  two  external  peripheral 
devices.  The  data  is  sent  to  a  Data  Collection  Platform 
Radio  Set  buffer  for  eventual  transmission  over  the 
satellite  link  and  it  is  sent  out  to  a  monitor  for  visual 
display  of  the  sampled  signal  levels. 

The  primary  decision  in  the  development  of  the 
program  Digitize.com  is  the  duration  of  the  sampling 
interval  used.  There  are  two  factors  which  must  be 
satisfied  in  order  to  make  the  decision.  The  first  is  to 
ensure  that  enough  samples  are  taken  to  accurately 
represent  the  signal  levels.  The  desired  outcome  is  to 
define  long  term  variations  in  the  background  noise 
levels  of  the  two  VHF  frequencies  of  interest,  therefore 
the  sampling  may  be  spaced  at  intervals  far  enough  apart 
to  provide  an  overall  trend,  not  necessarily  every 
instantaneous  fluctuation.  The  second  factor  involved  in 
determining  the  sampling  interval  is  available  data 
storage  space.  Since  there  are  256  quantization  steps  for 
an  8-bit  system,  the  analog  signals  are  digitized  and 
then  represented  at  each  sampling  by  a  three  digit  number 
(100  to  355).  For  convenience  in  processing  the  data  at 
a  later  time  the  data  is  stored  in  ASCII  format, 
requiring  3  bytes  of  information  for  each  sample  of  each 
signal  level.  With  ten  channels  to  be  sampled  at  each 
interval  this  becomes: 
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3  bytes  x  10  samples  =  30  bytes/sample 


The  limiting  factor  is  the  storage  space  available  on  a 
magnetic  floppy  disk.  A  double-sided  double-density 
floppy  has  80  tracks  on  it  for  use.  There  are  9  sectors 
per  track  and  512  bytes  per  sector.  The  first  sector  of 
any  disk  cannot  be  used  because  it  contains  formatting 
information.  Solving  for  the  available  storage  space: 

80  trks  x  9  sec/trk  x  512  bytes/sec  -  512  =  368,128  bytes 

If  a  two  byte  time  stamp  is  added  to  the  end  of  each 
sample  to  serve  as  a  sampling  period  spacer  there  will  be 
32  bytes  used  at  each  interval.  In  this  way  the  sectors 
will  be  evenly  filled  with  16  samples: 

30  bytes/sample  +  2  bytes/sample  =  32  bytes/sample 

512  bytes/sector  /  32  bytes/sample  =  16  samples/sector 

Solving  for  the  number  of  sectors  avafable  on  a  disk: 

368,128  bytes  /  512  bytes/sector  =  719  sectors 

The  number  of  32  byte  samples  which  can  be  taken  is: 
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368,128  bytes  /  32  bytes/sample  =  11504  samples 


or,  alternatively  solving  for  this  quantity: 

719  sectors  x  16  samples/sec tor  =  11504  samples 

The  final  decision  to  have  a  five  minute  sampling 
interval  is  a  compromise  between  the  amount  of  available 
data  storage  space  and  the  desired  degree  of  accuracy. 
This  interval  accurately  represents  the  expected  overall 
tendencies  of  the  background  noise  levels  and  allows  the 
following  quantity  of  time  for  data  storage: 

5  min/sample  x  hr/60  min  x  day/24  hrs  s  288  samples/day 

11504  samples  /  288  samples/day  =  39.94  days  of  space 

A  special  code  to  store  the  data  in  a  specific 
location  within  a  specific  sector  based  on  the  actual 
date  and  time  of  the  sampling  period  has  been  devised  to 
create  a  sequential  storage  method  for  the  data.  Because 
of  this,  the  time  limit  for  storing  the  data  on  a  disk  is 
limited  to  the  same  time  of  the  same  day  of  the  month 
following  the  insertion  of  the  data  disk  or  31  days, 
whichever  is  least. 

All  of  the  sectors  to  be  written  to  on  the  disk  are 
accessed  sequentially.  To  assign  the  current  sector  and 
location  within  the  sector  to  which  an  absolute  disk 
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write  operation  is  completed,  the  following  calculations 
are  made: 

N  =  sample  number 

N  =  12  *  hour  +  minutes/5  +  date*288 
S  =  sector  number 
S  =  N  /  16 

P  =  data  pointer  to  a  specific  location  in  the  sector 
P  =  M0D(N,16)  *  32 

The  date,  as  well  as  the  time  in  hours  and  minutes 
are  obtained  from  the  clock  chip  and  the  calculations  are 
made  for  the  write  operation  following  each  sampling 
interval.  As  a  result,  there  is  no  overwriting  of  data 
for  a  month.  Only  up  to  a  month  of  data  for  the  ten 
sampled  signal  levels  can  be  stored  on  a  disk.  It  must 
then  be  removed  from  the  noise  recording  system  disk 
drive  and  a  new  blank  formatted  disk  (using  PCDOS  2.1  or 
later  format)  inserted  to  collect  another  month  of  data. 
It  is  important  to  note  that  the  32  byte  data  samples  are 
stored  on  the  disks  by  an  absolute  write.  The  data  is 
not  in  a  file  structure  and  therefore  cannot  be  observed 
on  a  computer  at  this  time.  It  is  not  until  the  next 
step  in  the  data  processing  is  accomplished  that  the  data 
can  be  observed  by  the  user. 

2.  Datacopv.com 

The  program  Datacopy.com  as  shown  in  Appendix  B 
is  the  second  stage  of  the  data  processing  shown 
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previously  in  Figure  4.15.  The  program  is  written  using 
a  public  domain  version  of  the  higher  level  language 
FORTH  and  is  stored  on  an  IBM  PCDOS  2.1  format  disk  as 
an  executable  file.  The  FORTH  language  was  found  to  be 
the  simplest  and  quickest  means  of  accomplishing  the 
desired  outcome.  Also  stored  on  the  same  disk  is  an 
empty  data  file  called  Acquired.dat  which  will  eventually 
store  all  of  the  data  acquired  on  the  data  disk  by  the  RF 
noise  recording  system.  The  purpose  of  the  program 
Datacopy.com  is  to  perform  an  absolute  read  operation  on 
the  data  disk  and  then  copy  the  data  into  the  DOS  file 
Acquired.dat.  Once  the  data  is  in  the  DOS  format  file  it 
is  readily  available  to  the  user. 

The  transfer  operation  must  be  accomplished  on 
an  IBM  Personal  Computer  with  two  disk  drives  booted  with 
a  DOS  of  version  2.1  or  later.  To  perform  the  data 
transfer  the  following  steps  must  be  accomplished: 

a.  Place  the  disk  with  the  two  files  Datacopy.com 
and  Acquired.dat  on  it  into  disk  drive  "A" 

of  the  computer. 

b.  Place  the  data  disk  retrieved  from  the  noise 
recording  system  woth  up  to  a  month  of  data 

on  it  into  disk  drive  "B"  of  the  IBM  computer. 

c.  Type  in  the  command  "Datacopy  Acquired.dat'  at 
the  system  prompt.  The  program  will  direct  the 
user  from  this  point. 

Data  is  transferred  from  the  data  disk  in  IK 
blocks  and  the  program  continuously  indicates  how  many 
"K"  of  data  have  been  transferred  into  the  file 
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Acquired.dat.  The  transfer  can  be  stopped  at  any  time  by 
pressing  any  key.  If  a  key  is  pressed  the  user  will  be 
left  in  the  FORTH  environment.  To  return  to  DOS,  type  in 
the  command  "bye".  To  transfer  a  complete  month  of  data 
(278K)  takes  about  12  minutes.  Once  the  operation  is 
completed  and  the  program  has  stopped  transferring,  press 
<CTRL>C  or  <CTRL>  ALT  DEL.  The  entire  contents  of  the 
file  Acquired.dat  can  be  observed  by  entering  the  command 
"type  Acquired.dat",  The  data  is  now  stored  in  a  DOS 
file  format  and  the  user  is  ready  to  proceed  with  the 
next  step  of  the  data  processing. 

3.  Micro  to  Mainframe  File  Transfer 

To  allow  for  greater  versatility  and  storage 
capability,  the  DOS  file  Acquired.dat  created  by 
Datacopy.com  is  sent  to  the  IBM  3033  host  computer  using 
the  3270  terminal  mode  file  transfer  technique.  The  3270 
terminal  mode  file  transfer  is  a  method  of  file  transfer 
that,  can  only  be  used  when  the  PC  is  operating  as  a  3276 
model  2  display  terminal  and  is  connected  to  an  IBM  host 
computer  running  either  SIM3278/VM  or  SIM3278/VTAM  [Ref. 
13:p.  8-73.  To  perform  a  3270  terminal  mode  file 

transfer  the  IBM  PC  to  be  used  must  be  connected  to 
either  SIM3278/VM  or  SIM3278/VTAM  on  the  host  computer. 
A  number  of  IBM  personal  computers  with  this  capability 
are  located  at  the  Naval  Postgraduate  School.  Reference 
13  provides  complete  documentation  on  the  file  transfer 
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process . 


Using  a  terminal  meeting  the  required 
specifications  the  command  "send  Acquired.dat  Acquired 
Data  A1  (ASCII  crlf".is  entered.  The  disk  with  the  file 
Acquired.dat  containing  the  signal  level  data  must  be  in 
the  PC  disk  drive.  The  PCDOS  file  Acquired.dat  is 
transferred  to  the  users  "A"  disk  into  a  file  designated 
as  FN  FT  Acquired  Data.  The  data  will  be  sent 
sequentially  to  the  file  in  the  order  of  the  days  of  the 
month,  i.e.,  from  the  date  01  to  the  date  31.  If  during 
any  of  these  days  there  was  no  data  collected,  the  file 
will  be  blank  at  that  point. 

4.  Convert  Fortran 

Upon  the  completion  of  the  microcomputer  to  host 
computer  file  transfer,  the  next  step  in  the  processing 
of  the  data  is  to  rearrange  or  "convert"  the  data  in  the 
newly  created  acquired  data  file  into  a  more  convenient 
format.  The  program  Convert  Fortran  as  shown  in  Appendix 
C  is  designed  to  accomplish  this.  Once  the  data  is 
rearranged  by  this  program  it  can  be  used  as  a  hardcopy 
sequential  record  of  the  signal  level  data.  The  new 
format  will  also  allow  for  plotting  of  the  data  using 
DISSPLA  (Display  Integrated  Software  System  and  Plotting 
Language) . 

One  month  of  data  transferred  to  the  main  frame 
file  Acquired  Data  is  approximately  3900  lines  in  length. 
The  30  character  segments  of  data  (10  channels  of  data 
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with  3  digit  characters  used  to  represent  each  signal 
level,  i.e.,  100  to  355)  are  written  consecutively  into 
the  80  character  per  line  file.  Each  30  character  data 
segment  is  separated  by  a  two  character  spacer  (00,  01, 
02,  03,  04,  or  05).  The  program  Convert  Fortran  will 
read  the  30  character  segments  one  at  a  time  out  of  the 
Acquired  Data  file,  assign  a  sequential  number,  then 
write  them  to  the  new  file.  The  following  steps  must  be 
followed  to  accomplish  the  rearrangement: 

a.  Change  the  FN  FT  of  Acquired  Data  into  FN  FT  File 
FT08F001  and  save  this  newly  named  file. 

b.  Ensure  that  the  variable  named  ISTART  in  the 
program  Convert  Fortran  is  assigned  the  correct 
numerical  value.  ISTART  is  determined  by  the 
amount  of  data  characters  from  the  last  30 
character  segment  in  the  second  row  of  the  data 
file  Acquired  Data  which  have  not  overflowed  with 
the  rest  of  the  segment  onto  the  next  line.  An 
example  is  presented  in  the  system  performance 
section  of  this  thesis. 

c.  Run  the  program  Convert  Fortran  and  save 
the  rearranged  datafile  by  issuing  the 
following  sequence  of  commands  at  the 
appropriate  times: 

-  "Fortvs  Convert" 

-  "Record  on" 

-  "Load  Convert( start" 

-  "Record  off" 

d.  Temporarily  store  this  new  sequentially  arranged 
data  in  any  FN  FT.  Delete  any  extraneous  entries 
from  the  file  so  that  only  the  sequential 
numbering  and  data  remains  in  it. 

e.  Erase  the  file  previously  designated  File  FT08F001 
in  Step  a.  (the  original  is  still  in  the  file 
Acquired  Data). 
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f.  Change  the  FN  FT^of  the  file  created  in 
Step  d.  to  File  FT08F001 . 

Once  these  steps  are  completed  the  user  has  a 

sequentially  arranged  file  of  the  30  character  segments 

obtained  during  each  sampling  interval.  A  copy  of  this 

can  be  obtained  by  the  command  ’•print".  The  data  in  File 

FTOOFOOl  is  also  prepared  for  the  final  stage  of  the  data 

processing . 

5.  Earfaraf  Fortran 

In  this  final’  step  of  the  data  processing  a 
graphical  representation  of  the  data  in  File  FT08F001 
created  by  Convert  Fortran  is  obtained.  The  user  must  be 
located  at  a  terminal  with  an  attached  Tektronic  618 
graphics  monitor  to  run  the  program  Eqrfgraf  Fortran 
shown  in  Appendix  D. 

Once  the  necessary  minor  modifications  are  made  to 
the  program  Egrfgraf  Fortran  to  provide  tne  desired 
output  ,  the  following  commands  are  issued: 

"Fortvs  Eqrfgraf" 

"Disspla  Eqrfgraf" 

A  plot  of  the  selected  channel(s)  cf  data  is  displayed  on 
the  Tektronic  6 1 8  monitor.  A  hard  copy  of  the  graphs  may 
be  obtained  at  this  time. 

D.  SYSTEM  PERFORMANCE 

A  listing  of  selected  calibration  settings  e 

input  radio  frequency  signal  strengths  in  microvolts  to 

the  corresponding  analog  to  digital  converter  input 
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voltage  levels  and  quantized  signal  levels  is  shown  in 
Table  V.  The  system  was  placed  into  operation  in 
November  of  1986.  The  input  signals. were  sampled  and 
recorded  on  magnetic  disk.  The  data  recorded  directly 
onto  the  disk  was  transferred  to  the  DOS  file 
Acquired.dat.  Completing  the  next  required  step  of  data 
processing,  the  DOS  file  was  transferred  to  the  main 
frame  and  stored  in  the  file  Acquired  Data.  A  copy  of  a 
portion  of  this  file  is  shown  in  Figure  4.16.  The 
program  Convert  Fortran  rearranged  the  data  into  the 
sequential  format  file  callea  File  FT08F001  as  shown  in 
Figure  4.17.  The  final  step  of  data  processing  was 
completed  by  running  the  program  Eqrfgraf  Fortran  and  the 
plots  obtained  are  shown  in  Figures  4.18  to  4.27. 


78 


TABLE  V 


SYSTEM ' PERFORMANCE  LEVELS 


Input  RF  signal 
level  to  the 
receivers 
( mlc  royoi t5 ) 

Corresponding 
input  voltage 
level  to  the 

A/D  converter 
(volts) 

Corresponding 
quantized  signal 
level  stored 
in  data  file 
(unitless) 

0.0 

0.00 

100 

1.0 

0.32 

116 

2.0 

1 .00 

151 

3.0 

1  .77 

190 

4 .0 

2.49 

227 

5.0 

3.07 

257 

6.0 

3.57 

282 

7.0 

3.98 

303 

8.0 

4.36 

322 

9.0 

4.68 

339 

10.0 

5.00 

355 
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3"  JiV-’  j!  I"  1  *a  1  =  1  105!  5=2.751552*' o 
56  22  ’27j*3.'1  >>0:1  1’51  7624635  1271 
=  7  2  26  2  ’  J  ’.0*)1  =  3  ' =  3  l1  61  c  = 

5  3  11  =7.7  =?21  l  '-c  "5  1  l  J31-;  1240  05  =27J 

50  7  3  >-772^31  =67.33  l  7  6 1  s  =  0  3o05*it  i 
ol  2  2*0  751  4'-  iibOil  I3*l6s25i25625c 
61  22  0  pilJM  >C.-  5*  10  r  1  >1  * -.*♦  2  3?  2"  3 

6  2  223  2:07  30  1  1j755  t'OOl  70  2-1  05  427  1 

6"  >2  =27  >2311  =->  ’54  l  2  0  1  '*  0  2  4  5  2  5  3  "2  7  0 


Figure  4.17  Sample  File  FT08F001  File 
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CHANNEL  01  SIGNAL  LEVELS 
FOR  14  -  17  NOV  1986 
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CHANNEL  03  SIGNAL  LEV 
FOR  M  -  17  NOV  198 


Figure  4.20  Channel  03  Signal  Levels 


CHANNEL  04  SIGNAL  LEVELS 
FOR  H  -  17  NOV  1986 
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CHANNEL  05  SIGNAL  LEVELS 
FOR  M  -  17  NOV  1986 


Figure  4.22  Channel  05  Signal  Levels 
86 


SAMPLE  NUMBER/TIME 


CHANNEL  06  SIGNAL  LEVELS 
FOR  14  -  17  NOV  1986 
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CHANNEL  07  SIGNAL  LEVELS 


Figure  4.24  Channel  07  Signal  Levels 
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CHANNEL  08  SIGNAL  LEVELS 
FOR  14  -  17  NOV  1986 


Figure  4.25  Channel  08  Signal  Levels 
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CHANNEL  09  SIGNAL  LEVELS 
FOR  14  -  17  NOV  1986 


Figure  4.26  Channel  09  Signal  Levels 
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CHANNEL  10  SIGNAL  LEVELS 
FOR  H  -  17  NOV  1986 


Figure  4.27  Channel  10  Signal  Levels 
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V.  RECOMMENDATIONS  AMD  CONCLUSIONS 


A.  RECOMMENDATIONS 

There  are  several  avenues  of  research  to  be  pursued 
which  could  provide  further  correlation  between 
precursory  radio  frequency  emissions  and  earthquaKe 
activity.  By  utilizing  the  present  design  of  the  VHF 
noise  measurement  and  recording  system,  several 
modifications  to  the  hardware  and  software  may  provide 
more  explicit  results.  Limitations  would  of  course  be 
related  to  the  available  funding. 

The  range  of  frequencies  to  be  monitored  for  activity 
could  be  vastly  increased  by  the  introduction  of  a 
scanner  receiver  in  place  of  the  ten  fixed-range 
receivers.  There  are  a  large  number  of  commercial  grade 
scanner  receivers  capable  of  measuring  signals  in  the  KHz 
range  to  as  high  as  the  GHz  range.  Many  of  then,  are 
fully  programmable,  allowing  the  user  to  select  a  large 
number  of  channels  to  be  scanned.  The  scanning  of  the 
signals  at  the  selected  frequencies  would  need  to  be 
synchronized  to  the  sampling  rate  of  the  analog  to 
digital  converter.  This  would  require  the  design  of 
external  digital  circuitry.  By  incorporating  a  log- 
periodic  antenna  array  and  a  wideband  amplifier  into  the 
system  with  the  scanner  receiver,  a  much  wider  band  could 
be  monitored  for  radio  frequency  activity. 
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An  increase  in  the  measurement  accuracy  of  the  VHF 
signal  levels  being  monitored  could  easily  be  achieved  by 
increasing  the  sampling  rate.  A  simple  software 
modification  to  the  program  Digitize.com  would  accomplish 
this.  However,  the  increased  sampling  rate  would  also 
require  having  a  much  larger  data  storage  capability. 
The  use  of  hard  disk  drives,  or  even  magnetic  bubble 
memory  units  for  data  storage  would  accommodate  the 
increased  signal  sampling  rate,  but  the  expense  involved 
must  be  considered. 

Now  that  the  recorded  data  is  in  a  very  manageable 
form,  the  potential  for  much  more  detailed  long  range 
data  analysis  exists.  Programs  to  examine  the  recorded 
signal  level  data  could  be  created  to  perform  a  number  of 
mathematical  operations.  It  will  be  important  to  search 
the  data  over  a  long  period  of  time  to  deternine  routine 
patterns  and  cycles.  Once  a  firm  data  base  has  been 
established,  the  ability  to  distinguish  any  unusual 
disturbances  may  provide  a  deeper  understanding  of  the 
phenomena  being  observed. 

D.  CONCLUSIONS 

The  design,  construction  and  implementation  of  the 
computer  controlled  radio  frequency  noise  measurement  and 
recording  system  for  detection  in  the  30.45  KHz  and 
150.75  MHz  range  provides  a  means  of  analyzing  acquired 
data  rapidly  and  accurately.  The  interfacing  of  this 
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system  to  a  digital  satellite  link  adds  a  reliable  remote 
sensing  capability.  Future  analysis  of  the  data  may 
enhance  the  potential  for  using  radio  frequency  emissions 
as  a  precursor  to  future  earthquake  activity. 

At  the  present  time  more  research  is  required  to 
establish  a  correlation  between  radio  frequency  emissions 
and  earthquake  activity.  To  overcome  the  credibility  gap 
which  presently  exists  in  the  U3e  of  RF  emissions  as  a 
precursory  device,  a  detailed  documentation  of  the 
acquired  data  must  be  accomplished.  In  this  way,  a  basis 
for  future  predictive  capability  can  be  established . 
Several  criterion  need  to  be  identified  to  determine 
radio  frequency  precursory  anomalies.  First,  the  day  to 
day  normal  activity  pattern  must  be  clearly  and 
quantitatively  defined  by  extensive 

sampling  of  the  background  noise  activity.  Second,  the 
anomalous  RF  activity  has  to  be  defined  quantitatively 
arid  then  evidence  shown  that  the  "anomaly"  was 
significant  or  unique.  Third,  a  physical  or  casual  link 
between  the  RF  anomaly  detected  and  an  actual  earthquake 
must  be  established.  In  this  way  the  radio  frequency 
emissions  may  then  be  absolutely  defined  as  precursory. 
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MAIN  PROGRAM t 


1000  <303 UV  9000tREN  INITIALIZATION  ROUTINES 

1100  OOSUS  2000»REH  WAIT  FOR  3  MINUTE  DATA  TIHC 
1120  003U8  JOOOiRZN  0«T  10  1AMPLE3 

1140  603US  4000 iRCM  CONVERT  10  SAMPLES  TO  ASCII 

1200  RCADO.Sll  REN  READ  SECTOR  FROM  DISK  INTO  7*NORV 

1220  003 Ut  5000  i REN  PUT  NEW  ASCII  IN  BUFFER 

1S00  MRtTCCS.SliREM  WRITE  IUPFCR  1ACK  TO  DISK 

1SOO  003UI  7000 1 REN  VERIFICATION  OF  SAMPLING  CRITERION 

UO0  OOSUS  1000  1  REN  OUTPUT  CHANNEL  OATA  TO  TERNINAL/O.C.P  .R 

1900  OOTO  1 100 1  REM  LOOP 

2000  REN  -  WAIT  FOR  I  HI  MUTE  INTERVAL  AMO  CALCULATE  CURRENT  SECTOR 

2100  T ■CLOCKREAD i REN  OCT  MEN  TINE 

2120  IP(T-U)  OOTO  2100iREH  IS  OLO  TIME -MEM  TINE 

2140  U-T«REH  UFOATE  OLO  TINE 

2200  IF(MOQCT.S)OO)  OOTO  2100 tREM  NOT  S  MINUTE  INTERVAL 

2300  M*MCO(T.10O)«REN  CALCULATE  Mt MUTES 

2310  1PCM-40J  MTO  2100  PREVENTS  DOURLE  SAMPLE  ON  HR 

2320  N-CT-M)/100«AtM  CALCULATE  HOURS 

2400  0-(NP(*C20«)  ANO  IS  «  10>(  ir»(SS209 1  ANO  IS) » REN  OSTAIN  DATE 

2100  N»H«l2*M/S*0»:tSiREH  CALCULATE  SAHFLC  NUMtER 

2420  S-N/l 4 iRCM  SECTOR  MUHSCR 

2440  F -MOD  ( N •  1 4  ) ■  32 1 REM  POINTER  TO  OATA  IN  SECTOR 

2990  RETURN 

3000  REM  •  OCT  tO  SAMPLES  INTO  A  ARRAY 
SOSO  A-R1000 

3100  FOR  IO  TO  9 1  REM  10  SAMPLES 

3120  OUT (330. 1 )  1  REN  SET  UP  MUX  AMO  START  C»V 

3140  A< I )• INP(SSO) 1 REM  READ  RESULT 

3140  NEXT  I 

4000  REM  -  CONVERT  SAMPLES  TO  ASCII  PORHAT 

4100  FOR  1-0  TO  4l»RIM  FILL  1UPFER  WITH  42  SPACES 

4120  POKE(C»*20 )  1  REM  ASCII  SPACE 

4t40  NEXT  I 

4200  POR  1-0  TO  9 1  REM  NON  CONV  ANO  MOVE  10  SAMPLES 

4240  Q£N(C*I*4 » CAC  13*100 )  iREM  INTO  C  SUFFER  (RAMOf  *  100*311) 

4240  NfXT  I  sRCM  FINAL  <CR>  SflU.  IN  SUFFER 

4990  RETURN 

S000  REM  -  MOVE  MfW  SAMPLES  INTO  DISK  SUFFER  AREA  IM  RAM 

1101  FOR  !■!  TO  S 

1101  DPOKC(8*P*<  I»1 )  «DPEEK(CM  ) ) 

S1Q2  NEXT  I 

1103  FOR  l«S  TO  7 

1104  DPOKECB*PMI-2).OPCEK(CM>) 

SIOS  NEXT  1 

5104  FOR  1-9  TO  It 

1107  DFOK«(S*F*<I-3).0FfEKCC»I)l 

5105  NEXT  I 

1109  FOR  1*13  TO  IS 

SUO  DPOKCCl*P*<I-4),OPEEK(C*m 

Sill  NEXT  I 

SI  12  FOR  1-17  TO  19 

S1I3  OPOKE(S*P*(  I-§)  >OFEEX(C*I )) 

SI  14  NEXT  1 
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till  POO  1*21  TO  :s 

Silt  0PQXCC1*P*U-4).QPCEXCCM)) 

1117  NCXT  l 

Silt  FCX  l«2S  TO  27 

1119  DPOKCli*P*(l-7).DPCEX(C*in 

1129  NCXT  : 

S121  POO  |«29  TO  SI 

1122  DPOKf (l*P*t X«l).DPCIKCC*r>) 

1123  NCXT  l 

1124  POO  l»SS  TO  :s 

5125  OPOXCC«*P*<I-9).OPCCK(C‘I)) 

5 12 6  NCXT  I 

1127  POO  t«S7  TO  31 

1121  DPOKCCl*P*<I-10).DPCEX(C*m 

1129  NCXT  I 

1130  OCWI*P*30.M) 

1910  OCTUON 

7000  C  •  CLOCKOC/D 

7010  P  •  10»C  INP(*xC20S)  AND  IS)  •  tNPC*C20A)  AND  IS 

7020  0  •  l0»(tNP(*C2Q9)  AND  11)  •  INPCXC20*)  AND  IS 

7030  V  •  10MINP(*I200>  ANO  IS)  •  INP(U20C)  ANO  IS 

7100  POINT  *1  IIT  A/O  SIONAL  LCVCLS  POO « -  .P .  -  /“.0.“  /".V 

7200  OCTUON 

•000  POO  I»0  TO  9iOCH  5CN0  DATA  TO  MONITOO 

•010  POINT  "CHANNCL".!*!."  DIOITKU  SIGNAL  LCVCLi  AC  l  IMOO 
•Oil  POINT  -  - 
•020  NCXT  I 

•030  POINT  - . . . * . * . . . . 

•010  COMINtTC 1 10 »2«0 • 1 ) lOCN  INITIALIZI  OUTPUT  POOT 

•Oil  POO  l  ■  0  TO  9 

•017  Z  ■  AC  1 3  SIND  SINAOV  DATA  OUT  OP 

•  040  V  >  CQHOUTCZ)  tOCM  POOT  4  CP4)  TO  THC  O.C.P.O.J. 

•041  NCXT  I 
•100  OCTUON 

9000  OCM  -  INITIALIZATION  CODC 

9100  U'OiOCM  SCT  OLD  TINC  TO  IHPOSSIfLC  VALUC 

9120  C«<1700«OCN  ASCII  CON V CO SI ON  SUPPCO 

9140  ••CllOOiOCM  IUPPCO  POO  OISK  DATA 

9990  OCTUON 
S20Q0  STOP 
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APPENDIX  B 


DATACOPY.COM  COMPUTER  PROGRAM 


•  PROGRAM  TO  If  USED  ON  AN  IBM  MICRO  • 
»  TO  00  AN  ABSOLUT!  SECTOR  READ  PROM  « 

•  THf  MAONCTIC  DISK  AND  THEN  STORf  • 

•  TMt  OATA  IN  A  DOS  PORMAT  PILE.  ■ 


/  DATACOPV  PROORAM 

l  J  ‘THRU  /  LOAD  IN  THC  PROORAM  CODE 

•  DO- IT  IS  BOOT  /  SET  UP  POP  AUTO-RUN 

SAVE- SYSTEM  BiOATACORY.COM 

/  DISX  INTERPACE  VIA  BtOS 

CODE  INTIS  (S  OFFSET  SECTOR*  OP  —  P  J 

01  POP  AX  POP  OX  OX  XOR  720  ■  IX  MOV  BX  OIV  AX  PUSH 

QX  AX  MOV  OX  OX  XOR  BX  SNR  BX  DIV  BX  POP  AC  BH  MOV 

OX  AX  MOV  I  •  CL  MOV  CL  OIV  AM  AL  XCH9  AX  INC 
AX  CX  MOV  IX  OX  MOV  BX  POP  OS  PUSH  IS  POP  01  AX  MOV 
AM  AC  XCMO  1  •  AL  MOV  SI  PUSH  BP  PUSH  It  INT 
9  •  AX  MOV  U<  IP  AX  OEC  THEN  BP  POP  SI  POP  IPUSH 
ENO  COOC 

CODE  SETUP  (S  BUFFER  OP  —  0FF1  SEC1  OP  OFFS  SECT  OP  ) 

OX  POP  IX  POP  0  (IX I  AX  MOV  4  (IX I  CX  MOV  AX  SMC 

CX  PUSH  AX  PUSH  DX  PUSH  S 12  •  CX  ADD  AX  INC 
CX  PUSH  AX  PUSH  OX  PUSH  NEXT  ENO  CODE 
«  RBIOS  2  SETUP  INTIS  >R  INTIS  R>  OR  OISK-CRROR  !  I 

«  MHOS  S  SETUP  INTIS  >R  INTIS  R>  OR  OISK-ERROR  f  I 

i  BIOS- 10  l  *  J  RBIOS  IS  RCA  0- BLOCK  Cl  MB  I  OS  IS  WRITE- BLOCK  l 

i  SET- 10  /  SET  UP  TO  READ  RAN.  WRITE  TO  DOS 

(  DOS  |  PILE- 10 

I  FORTH  |  CM  RBIOS  IS  REAO-BLOCX  » 

»  DO-COPY  /  COPY  SI  OAYS  WORTH  OP  OATA 

2?t  0  DO 

I  S4f  •  /  DATA  STARTS  ON  SECTOR  IB  OP  OATA  OISX 

IS  CHIT  t  I  .1  /  (TWO  SCCTQRS/K.  IK/ SCREEN  ) 

I  COPY  OISK-ERROR  I 
KEY?  OR  IP  OUIT  THEN  LOOP  I 

i  PROMPT  /  WAIT  FOR  THE  USER  TO  BE  READY 

."  INSERT  OATA  DISK  IN  ORIVE  It.  DOS  OtSK  IN  ORIVE  A.  -  CR 
PRESS  ANY  KEY  WHEN  READY.  M  CR  KEY  DROP  t 

l  DO- IT 

CR  PROMPT  /  GET  THE  DISKS  IN  PROPERLY 

DEFAULT  /  OPEN  THE  OUTPUT  PILE 

SET- 10  /  SET  UP  TO  READ  SATA.  WRITE  TO  DOS  PILE 

DO-COPY  /  TRANSFER  THE  OATA  TO  THE  PILE 

CR  TRANSFER  COMPLETT" 

BYE  I  /  ANO  EXIT 
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APPENDIX  C 


CONVERT  FORTRAN  COMPUTER  PROGRAM 


C  . . . . . . . 

C  ■  THIS  PROGRAM  MIU.  CONVERT  *NC  DATA  FILE  THAT  MAS  TRANSFERRM  * 
C  «  FA  OH  AN  IIM  AC  TO  THC  MAIN  FRAME  INTO  A  NEW  FORMAT .  THIS  ■ 
C  *  NCM  FORMAT  HILL  ALLOW  FOR  FLOTTINO  OF  THC  DATA  USING  DISSPLA- ■ 


THC  DATA  ARRIVES  ON  THC  MAIN  FRAME  IN  THC  FILI  'ACQUIRED  DATA* 

FA  OH  THC  FCDOS  FILI  CALLID  *  ACOUIRM.DAT’  AS  APPROXIMATELY  J*00 
LIHCJ  OF  OATA.  THC  SO  CHARACTER  SEGMENTS  OF  DATA  (10  CHANNELS  OF 
DATA  WITH  S  CHARACTERS  Ft*  CHANNCU  ARE  WRITTEN  INTO  THC  SO 
CHARACTER  FIX  LIHC  PILE.  SEPARATED  It  TWO  DIOIT  T1HC  STAMFS  U.C. 
01  *02 . OS . 04 .  OR  OS).  THIS  PROGRAM  FIAOS  THC  SO  CHARACTER 
SCOHCNTS  ONC  tV  ONC  OUT  OF  THC  ACOUIRED  DATA  PILE  AND  ARRANGES 
THCM  IN  A  SCOUCNTIAL  FILI.  IN  THIS  WAV  THC  DATA  IS  AVAILABLC 
IN  AN  EASILY  OISCCXNAILK  FOAM  FOX  THC  USER  ANO  ALSO  IS  RCAOILY 
GRAPHED.  MINOR  AOJUSTHCNTS  IN  THIS  PROGRAM  HAY  M  NECESSARY 
TO  ACCOUNT  FOX  VARIATIONS  IN  THC  MICRO  TO  NAIM  FRANC  TRANSFIX 
FXOCCSS.  THC  FOSS  I  ILK  CMANOCS  ARK  NOTtD  WITHIN  THC  FROORAM. 
FOLLOW  THC  STIRS  OUTLINCO  KLOW  TO  ARRANQC  THC  OATA  IN  A  FORMAT 
FOR  USC  IN  A  FROORAM  TO  OR  ARM  I  CALL  Y  OISFLAY  THC  OATA. 


STCF  1 1  FUT  THC  DATA  IN  THC  FILI  ’ACOUIRED  DATA*  INTO  A  FtLC 
CALLID  'FILI  FT0CF00I*  AFTCR  ANY  NON  ESSENTIAL 
INFORMATION  THAT  MAS  TRANSFORM  HAS  SCIN  MITM  OUT. 
STIF  2 1  MAXI  ANY  NECESSARY  CHANCES  IN  THIS  FROGRAM.  ISFECIALLY 
TO  THC  “1ST ART"  FARAMCTO. 

STEP  Si  RUM  THIS  FROORAM  ANO  SAVE  THC  NCNLY  FORMATTED  OATA 

FILE  «Y  ISSUING  THC  FOLLOWING  SEQUENCE  OF  COMMANDS  AT 
THC  AFFROFRIATC  TIMES i 
'FOR TVS  CONVERT’ 

•RECORD  ON’ 

•LOAO  CONVERT (START* 

•RECORD  OFF’ 

TEMPORARILY  STORE  THIS  DATA  IN  ANY  FILENAMC.PILCTYRC. 
STEF  4 1  ERA.’!  THE  FILE  CALLED  ’FILE  FT0SF001’  CREATED  IN  STEF  1. 
STEF  Si  CHANGE  THC  NAME  OF  THE  FILE  CREATED  IN  STEF  4  INTO 
'FILE  PTOSPOOI * • 

STEF  «i  GO  TO  THE  FILE  ’EQRPGRAP  FORTRAN’  ANO  PROCEED. 
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IMT1QCR  1  .J.K.L. (START 
OUR  AC  TOM  A(UO).ICJO) 

C  MOTH  ta  INSTRUCTIONS  ON  NON  TO  SC7  THU  FARANfTIR  "(START**! 
I START  •  U 

L  •  1 

I  RSAOCS  • IS itNO  >  1000  )  A 
IS  FORMAT (S0A1/S0A1 ) 

J-0 

DO  20  I  •  (START* 1. SO 
20  Kit  ■  AC  I  •(START) 

WRITS (4*  SO)  L.l 
SO  FORMAT (T2 .  (S.1X.IIA1) 

U  •  L  *  I 
J  ■  St  -  (START 
00  SO  K  •  1.4 
DO  40  I  ■  1.10 
40  SCI)  •  A( (*J) 

MR(Tt(4.S8)  L.S 
L  •  L  •  1 
SO  J  ■  J  •  12 

00  40  I  •  (.(START 
40  1(1)  •  A( 140  •  (START  •  () 

00  TO  S 
1000  CONTI HUS 
A STURM 
KNO 
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APPENDIX  D 


EQRFGRAF  FORTRAN  COMPUTER  PROGRAM 


c  . . . 

C  •  THU  PROGRAM  HILL  GRAPH  THf  DATA  STORED  IN  A  PILE  CALLED  'FtLI  ■ 
C  »  PTOiPflll '  CREATED  IV  THC  PROGRAM  ‘CONVERT  PORTRAN* .  YOU  MUST  »C  • 
C  *  AT  A  TERMINAL  WITH  AN  ATTACHED  TEXTRON  I C  All  TO  RUN  IT. 

C  . . . . . . . . . 

c 

c 

c 

C  ISSUE  THE  FOLLOW  I  NO  COMMANDO  FROM  CMS  TO  RUN  THIS  PROGRAM  I 
C  'FOR TVS  IQRPORAF’ 

C  ' 0IS2FLA  ICRFGRAF* 

C  A  GRAPH  OF  THE  SELECTED  COLUMNS  OF  OATA  WHICH  HAVE  NOT  SEEN 
C  COMMENTED  OUT  IELOU  WILL  BE  DISPLAYED  ON  THE  TEXTRON  I C  All 
C  MONITOR.  HARO  COFY  MAY  1C  OBTAINED  AT  THIS  TIME. 

C 

C 

. . . . . . ••• . . . 

C 

c 

c 

INTEGER  I 

C  MAKE  SURE  THESE  ARRAYS  ARE  LARGE  ENOUGH  TO  STORE  THE  DATA! 

REAL  TC020). 00(020), 01 (020). 02(129). DJI020).  04(020) 

REAL  09(1201  •  DA  (ISO)  *07(120)  ■  00  (ICO)  .09(120 1 
C  SET  THC  INOCX  (I)  ON  THIS  DO  LOOP  TO  THE  CORRECT  RUMKRl 

C  I.E.*  TO  THC  NUMBER  OF  ROWS  (SAMFLES)  TO  BE  READ  FROM  THE 

C  DATA  FILS i  (NOTES  THERE  ARE  2tl  SAMPLE S/CMANNCL/OAY ) 

DO  20  I'1.7«S 

REAO(0'A01)T(!).00(!).01<!).D2(n.DS(!).D4m.09(n.DA(n 
1.07(1). 00(1). 09(1) 

20  CONTINUE 

AO l  PQRMAT(ES.0< IX. 10(11.07) 

CALL  TEXAli 

CALL  FAOE  (U.O.t.l) 

CA<i.  0L0WUFU.2) 

CAU.  ARSA20  (f.I.A.O) 

CALL  YNAMC  ('SIGNAL  LEVEL I*. 100) 

CALL  XNAMC  ('SAMPLE  NUMVCR/TIMCS * .  100 1 
C  ENTER  THE  CHANNEL  MUMICR(S)  TO  BC  GRAPHED: 

CALL  HEAOIN  ('CHANNEL  91  SIGNAL  LEVELS* *  .100.2.2) 

C  ENTER  THE  CORRECT  DATEt 

CALL  HEAOIN  CPOR  14  -  17  NOV  191  At*  .100.S.2) 
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«  nwn  X  AM  V  RAMGCS  OP  THC  OPAPHi 
IP  (S.90.IIS.100.2Q.J40] 

10(1. i) 

.INC 

«  ENOINO  DATA  PACAMCTCR  ON  THf  OWKUO  SCLtCTVOi 

II 

'VC  (T.D0.7IS.I) 

Ci 

'VC  (T.P1.7IS.0) 

Si 

'VC  (T.O3.7CS.0) 
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